Characteristics of reaction-diffusion on scale-free networks 
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In this work we examine some characteristic properties of reaction-difTusion processes of the 
A + A ^ Q type on scale-free networks. Due to the inhomogeneity of the structure of the substrate, 
as compared to usual lattices, we focus on the characteristics of the nodes where the annihilations 
occur. We show that at early times the majority of these events take place on low-connectivity 
nodes, while as time advances the process moves towards the high-connectivity nodes, the so-called 
hubs. This pattern remarkably accelerates the annihilation of the particles, and it is in agreement 
with earlier predictions that the rates of reaction-diffusion processes on scale-free networks are much 
faster than the equivalent ones on lattice systems. 
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I. INTRODUCTION 

The A + A ^ Q reaction is one of the most fundamen- 
tal annihilation processes jwith many applications 1]. In 
a number of papers 0, Q] this reaction was studied on 
scale-free networks, both computationally and analyti- 
cally. In a recent Letter 0] we reported some unusual 
properties of reaction-diffusion processes taking place on 
scale-free networks. These networks represent a special 
class of random networks, where the probability that a 
node is connected to k other nodes of the network follows 
a power-law form: 



P{k) - k-^ 



(1) 



where the value of the exponent 7 usually lies in the range 
2 < 7 < 4. A large number of systems from markedly 
different disciplines have been found to fall in this class, 
and new systems are added continuously in the list 0, Q . 
This ubiquity explains the intense interest devoted to the 
study of the complex networks field. The unusual struc- 
ture of the scale- free networks has been shown to strongly 
modify the properties of dynamic processes, where par- 
ticles move from node to node along the existing links, 
as compared to the same heavily-studied processes on 
lattice systems or continuous space. 

Reaction-diffusion mechanisms are a frequently en- 
countered mechanism for kinetics, where particles diffuse 
in a space and react in a predefined way upon encounter- 
ing other particles. For example, during the A + A ^ Q 
process, particles of the A type diffuse and annihilate 
when they collide with other particles of the same type. 
In Euclidean lattices, the simplest mean- field approach 
predicts a linear increase of the inverse particle concen- 
tration pit) with time t. A large number of studies 
0i Hi Q J though, demonstrated that the density actually 
scales as 



P{t) pa 



(2) 



where / < 1, and po is the initial particle concentration 
at time t = 0. In a space with dimensionality d this 



exponent equals / = d/dc, for d < d^ and / = 1 for 
d > dc- The critical dimension for the A + A reaction is 
dc = 2. This 'anomalous' behavior has also been observed 
when the substrate of diffusion has different geometry, 
such as a fractal structure where / — dg/dc, and ds 
is the spectral dimension. 

The maximum value of / in all the above cases has been 
/ = 1. In scale- free networks, though, we recently found 
numerically that the process follows a different mecha- 
nism. For networks with 7 < 3.5, the concentration de- 
cay still follows a power law, which for a short interval is 
described initially by / < 1, but soon after that exhibits 
a crossover to a power law with an exponent / > 1. 
This is a clear indication for the rapid acceleration of the 
process, which we attributed mainly to the existence of 
the hubs. Following this, Catanzaro et al. Q developed 
a theoretical framework for 'uncorrelated' networks (i.e. 
networks with no degree-degree correlations, 
They predict a behavior 



see 
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2<7<3 
ln(t), 7 = 3 
t, 7>3 



(3) 



A detailed comparison between the two models was re- 
cently presented in Q , where it also became evident that 
the exact behavior of the particle concentration depends 
on details of the structure, such as, for example, the min- 
imum number of links on a node. 

In the present work we study the reaction-diffusion 
A + A ^ model on scale-free networks, with empha- 
sis on the nodes where the annihilation events take place 
on, and we try to better understand the mechanism that 
leads to this new type of behavior. 



II. THE MODEL 

The networks that we use in our simulations _are pre- 
pared with the standard configuration model 12]. Ini- 
tially, an integer k value representing the number of links 
is assigned from a random distribution obeying Eq. [T] to 
every one of the iV nodes in the network. Random pairs 
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of links are chosen and connect two nodes, but double 
links and self-connections are not allowed. This network 
creation method may create a certain number of sepa- 
rated clusters (depending on the value of 7). We isolate 
and use only the largest cluster of the network, where the 
particles are allowed to diffuse. 

During the reaction-diffusion process the number of 
particles M{t) in the system is reduced with time, and 
we denote with p{t) the particle concentration at time 
t. The initial particle concentration is po- A number 
of A^o — Po^ particles are placed on randomly selected 
nodes (in simulations of this work we use po = 0.5). 
Then, we randomly choose one particle and pick one of 
the neighboring nodes where this particle is located at. If 
this new node is empty, the particle moves and occupies 
its new position. If this node is already occupied, then 
the two particles annihilate and are immediately removed 
from the system. Time is advanced inversely proportion- 
ally to the current number of particles, by 1/M(t). We 
repeat this procedure by continuously selecting, moving 
and (possibly) annihilating particles. We perform 100 
different realizations of the walk on different networks of 
size 10^ nodes each. The largest cluster on these net- 
works is, of course, smaller and spans 35% to 100% of 
the system nodes, depending on the value of 7. 



III. RESULTS 

In order to gain insight into the reaction process, and 
given that the substrate of the diffusion is very inhomo- 
geneous, we try to understand the nature of the nodes 
where the annihilations take place. The connectivity of 
the nodes in the network varies largely from node to node, 
but the great majority has a small number of connections, 
which is the main characteristic of a scale-free network. 
A very small number of nodes are well-connected (the 
hubs), and the existence of the hubs has been shown to 
be the main reason behind most of the unusual properties 
reported for the scale-free networks. 

In this context, we have performed simulations of the 
A + A reaction-diffusion problem on scale-free networks 
of various 7 exponents. After t — 100 steps, when the 
particle concentration has diminished to small values, we 
calculate the percentage of annihilation events i^ann(^) 
that took place on a node with k links. Figure [T] shows 
-Fann(^) SiS SL function of k. The raw data for -Farm exhibit 
a kink at values of k around 100 to 1000, above which 
the curve increases monotonically with increasing k. We 
present one such curve in Fig. [T]for the case of 7 = 2.0. 
This feature (which is also size-dependent) is a numerical 
artifact, due to the rarity of nodes with very large degree 
k, and is removed when the data are binned. The four 
curves in the figure are the result of logarithmic binning 
at degrees fc > 10. 

In a mean-field treatment of a network with no degree- 
degree correlation we expect that the number of events 
would be proportional to the total number of links lead- 




FIG. 1: Probability, Fann(fc), for an annihilation event to take 
place on a node with k links as a function of k, during a 
process of t = 100 steps. Open symbols are averages over 
100 scale-free networks with 7=2.0, 2.5, 3.0 and 3.5 (top to 
bottom, shown in the plot), oi N = 10^ nodes each. The data 
have been logarithmically binned for k > 10. Lines represent 
Eq. |4l with the corresponding exponent 1 — 7 for each case. 
We also present unbinned data (dots) for the case of 7 = 2.0. 



ing to a fc-node, i.e. 

i^an„(fc) - kP{k) ~ fcl-T . 



(4) 



The exponent 1 — 7 is verified in all four cases, as can 
be seen in Fig. [Tl for the intermediate to large k regime. 
The result of Eq. 2] and the presented curves show that 
even though most of the events take place on low-fc nodes, 
which comprise most of the network, the probability for 
annihilation at the hubs is significantly larger than their 
relative appearance in a network. The behavior seems 
also not to be infiuenced by any degree-degree correla- 
tion in the network. The probability that nodes with 
low connectivity are connected to a hub can in fact be 
larger than the one predicted for a completely uniform 
distribution [ll[, due to the avoidance of double links 
and self-connectivity. 

The importance of the hubs can be made explicit if 
we ask what is the average number of events R{k) on a 
single node with k links, i.e. 



Rik) 



-Fann(fc) 
P{k) 



(5) 



According to Eq. [¥]we expect a linear increase with the 
number of links fc, since i?(fc) ^ k^^^ /k~^ ^ k. This 
behavior is verified in Fig. [21 at least for large fc values. 
The slope of 7 = 2 seems to be slightly different than the 
other curves, but this can be attributed to the increasing 
number of hubs at this 7 value, where the network has 
a very small diameter. The curves in Fig. [2] span many 
decades on the y-axis manifesting a very different role 
of nodes with different degrees. The case of 7 = 2 is 
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FIG. 2; Average number of annihilation events per node, 
R{k), with k hnks as a function of k, during a process of 
t = 100 steps. Curves are averages over 100 scale-free net- 
works with 7=2.0, 2.5, 3.0 and 3.5, of = 10*' nodes each. 



also the borderline value that separates networks with 
finite degree distribution moments from networks with 
infinite moments. The diameter of networks with 7 < 2 
is extremely small and the hubs connect practically the 
entire network within a few steps. In such networks with 
7 < 2, which are difficult to simulate numerically, we 
expect that the annihilation process will be exponentially 
fast and that almost all the reaction activity will take 
place on the hubs or on their immediate neighbors. 

In a typical walk oft — 100 steps about 2 x 10^ events 
are observed. If we focus on a particular node with a 
small degree we see that there is much less than one an- 
nihilation event occuring on this node during the process. 
On the contrary, more than 10^ events occur at the hubs 
having k ~ 10^ links. We, thus, conclude from Figures 
[1] and [2] that a hub is a much more active element in a 
network than any other node, although the majority of 
events still occurs over the large number of low-fc nodes. 

The results until now concerned the entire process af- 
ter averaging over the first t — 100 steps. We now turn to 
the time evolution of the diffusion-reaction mechanism. 
In Fig. [3^ we can see how -F'ann(fc = 1) varies with time. 
For all 7 values the probability starts initially with a cer- 
tain value and decreases monotonically with time before 
stabilizing to a much lower value. For 7 = 8 and 7 = 3.5, 
where hubs are not very strongly connected, the probabil- 
ity goes down to a value of 0.02 even for the longer times 
displayed. For 7 = 2 and 7 = 2.5, though, the proba- 
bility that an annihilation occurs on a node with fc = 1 
practically vanishes. This picture is reversed when we 
monitor the probability -Fknnlfcmax), i-e. the fraction of 
events taking place on the single most connected node of 
the network during the time interval [i, t+\] over the total 
number of events that occured during the same interval. 
Figure [5)3 suggests that now -F'ann(^max) starts with low 
values but continuously increases as time evolves. Note 




time 

FIG. 3: Time evolution of the average probability (a) 
fann(fc = 1) and (b) -F'ann(fc = fcmax), for uetworks with vary- 
ing 7 exponents. 

that for networks with 7 < 3 more than 10% of the events 
at long times take place on this node alone. 

In Fig. |4] we present the time evolution of the entire 
-Fann(^) distribution for different 7 values. The distribu- 
tion at i = 1 has a form similar to that of Fig. [T] When 
7 = 2, and as time evolves, the low-fc part decreases and 
the large-fc part increases. Asymptotically, the distribu- 
tion tends to uniform and for t ~ 40 there is no event 
taking place at nodes of small connectivity. For 7 = 2.5 
the picture is similar, but at i = 100 the distribution has 
not yet spread uniformly over the entire k range. For 
7 > 3 the right wing increases with a slower rate than 
when 7 < 3, and always remains significantly lower than 
the small- fc part, at least for the early time scale that we 
observe in our current work. 

Notice also that at time t = 100 the process on net- 
works with different 7 lies in different stages of the p{t) 
evolution, as can be seen in Fig. 1 of Ref. [2]. For 7 = 2 
the concentration has almost diminished to zero, while 
for 7 = 2.5 we are well inside the power law regime with 
/ > 1. When 7 = 3.0 we are located close to the crossover 
point between the two power law regimes, and for 7 = 3.5 
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FIG. 4: Time evolution of the probability distribution Fann(fc) 
for four different 7 values (shown on the plot) . Data have been 
logarithmically binned for k > 10. Each curve corresponds to 
different times displayed in the legends. 



the concentration has not reached the transition point 
and still decays slower than linear. 

In Ref. [3] we had shown that during an A + A ^ re- 
action particles tend to segregate. This is a very unusual 
type of behavior, since in regular Euclidean space the 
opposite effect (formation of depletion zone) is observed. 
We, thus, consider in more detail the cluster structure 
of the particles in the network. We have already seen 
(Fig. 3a of 0) that the fraction Qaa, defined as the 
number of the close contacts in the system over the to- 
tal possible number of contacts, increases as a function 
of time (except for 7 > 3), a clear indication towards 
particle segregation. In order to fully characterize the 
motion and the spatial distribution of the particles we 
now use the two-particle correlation function g{L), which 
we define as the number of particles at distance L over 
the number of particles at distance L when we assume 
an equal but uniform particle distribution. On random 
networks, this random distribution cannot be accurately 
estimated theoretically, as e.g. for simple lattices, due to 
the widely different local environment. In order to over- 
come this problem, for every network that we used and 
for every concentration encountered we placed a corre- 
sponding number of particles uniformly on this network 
and averaged over 20 different realizations. By calculat- 
ing the average number of particles at distances L we 
were able to determine the value of the denominator. 
Thus, if the particles are randomly distributed in the 
network this quantity will be 1. If the particles tend to 
occupy nodes far from each other this fraction will be 



FIG. 5: Two-particle correlation function g{L) as a function of 
the distance L for different 7 values (shown in the plot). The 
dashed line represents the uniform distribution limit g{L) — 1. 
Thicker curves correspond to earlier times. Times used in each 
plot: (a) t =3, 8, 20, 30. (b) t =6, 25, 100, 150, (c) t =12, 72, 
500, 850, and (d) t =15, 140, 1800, 3400. 



g{L) < 1, while g{L) > 1 shows that particles are clus- 
tered in close contact to each other. 

The first three plots in Fig. O (corresponding to 7 < 3) 
describe a picture where particles break their initial uni- 
form distribution and cluster with an increasing rate. 
The average number of neighbors within a few steps is 
significantly high and increases with time, while the num- 
ber of particles at longer distances are close to the mean- 
field assumption g{L) = 1. When 7 — 3.5 the correlation 
function is more similar to a typical A -I- A reaction on 
a normal lattice. Initially, a depletion zone is created 
around the A particles. This depletion zone is retained 
for the duration of the process, with the difference that at 
longer times we observe an increase of particles concen- 
tration at moderate distances. This increase is, though, 
much smaller than for the case of 7 < 3. 



IV. SUMMARY 

Summarizing, in the present study we focused on the 
details of the diffusion, such as the spatial organization of 
the particles and the time evolution leading to an annihi- 
lation reaction, with a scale-free network as the substrate. 

The majority of the particles annihilate on the low con- 
nectivity nodes. This happens mainly at early times, 
where the static distribution of particles on the network 
is the dominant factor. The particles annihilate in the 
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vicinity of their initial neighborhood, which on average 
includes mostly low-degree nodes (since the number of 
such nodes in the system is large). 

At later times, when a small number of particles re- 
main in the system, diffusion starts to play an increas- 



ingly important role. Particles encounter each other 
mainly on nodes where diffusion directs them to, i.e. on 
the hubs[l^. Thus, although the percentage of the hubs 
in the network is very small most reactions take place on 
them and they dominate the process at longer times. 
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